Transient Charging and Discharging of Spin-polarized Electrons in a Quantum Dot 
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We study spin-polarized transient transport in a quantum dot coupled to two ferromagnetic leads 
subjected to a rectangular bias voltage pulse. Fime-dependent spin-resolved currents, occupations, 
spin accumulation, and tunneling magnetoresistance (TMR) are calculated using both nonequilib- 
rium Green function and master equation techniques. Both parallel and antiparallel leads' mag- 
netization alignments are analyzed. Our main findings are: a dynamical spin accumulation that 
changes sign in time, a short-lived pulse of spin polarized current in the emitter lead (but not in the 
collector lead), and a dynamical TMR that develops negative values in the transient regime. We 
also observe that the intra-dot Coulomb interaction can enhance even further the negative values of 
the TMR. 



PACS numbers: PACS number 
I. INTRODUCTION 

A variety of new effects and novel devices have been 
reported during recent years in the context of the emerg- 
ing field of spintronicsJ^i 3 -^ One of the most challeng- 
ing milestones in this context is the development of a 
quantum computer, which would represent a great break- 
through in the processing time of certain mathemati- 
cal and physical problems £ In particular, the electron 
spin in quantum dots has been proposed building 
block for the implementation of quantum bits (qubits) 
for quantum computation^ 7 - An important recent de- 
velopment is the possibility to coherently control elec- 
tron states and electron spin in quantum dot systems 
with a precision up to a single-electron, thus demon- 
strating the feasibility of qubit implementation in a solid 
state system! 8 ! 9 ! 10 ! 11 ' 12 ! 13 Specifically, these experimen- 
tal realizations use high-speed voltage pulses to tunc 
the system levels in a coherent cycle for electronic ma- 
nipulation. Ac-driven quantum dot systems and dou- 
ble barrier structures have also been studied in the 
context of quantum pumps ; 14 ' 15 ' 16 ' 17 ' 18 supcrlatticesj^ 9 - 
Kondo effect ) 20 i 21 ' 22 i 23 and spin-polarized transport] 24 ' 25 
In addition to this, time-dependent transport has re- 
ceived growing attention in a variety of mesoscopic sys- 
tems that encompasses, to mention but a few, molecu- 
lar electronics , 26 ' 27 dissipative driven mesoscopic ring,— 
noisy qubits r 9 - and dynamical Franz-Kcldysh effect. 30 

In the context of spintronics a system of particular in- 
terest is composed of a quantum dot or a metallic is- 
land coupled via tunnel barriers to two ferromagnetic 
leads (FM-QD-FM). For example, in the nonequilibrium 
regime the following effects have been discussed: a spin- 
split Kondo resonance ; 31 ' 32 a spin-current diode effect^ 



zero-bias anomaly^ tunnel magnetoresistance (TMR) 
oscillations^ negative TMR^ spin accumulation^ and 
so on. In spite of all this activity, to the best of our 
knowledge, only very little work has been done on spin- 
polarized transport driven by ac-bias voltages! 24 ' 25 Here 
we study transient spin-resolved currents, occupations 
and TMR generated by a voltage pulse applied in one 
of the ferromagnetic leads. We use two complementary 
approaches to study the problem: nonequilibrium Green 
function (NEGF) and the Master Equation (ME). NEGF 
is used to give an exact solution in the noninteracting 
case, while the ME, valid in the limit fcgT ^> To (Fo 
is the characteristic level width), is used to demonstrate 
that the results obtained via NEGF are modified only 
quantitatively, not qualitatively, when Coulomb interac- 
tion is accounted in the sequential-tunneling limit. Both 
parallel (P) and antiparallel (AP) magnetization align- 
ments are considered. In the P case we find a magnitude 
and sign modulation of the spin accumulation in the dot, 
while in the AP alignment only the magnitude changes. 
For the current we observe a spike of spin polarized cur- 
rent in the emitter lead when the system operates in the 
P configuration. This effect gives rise to a dynamical 
negative-TMR just after the bias voltage is turned off. 



The paper is organized as follow. In Sec. II we de- 
scribe the formulation based on NEGF, and give explicit 
formulas for the noninteracting case. In Sec. III(a)-(c) 
we present numerical results based on Sec. II, and in Sec. 
111(d) we apply the master equation technique to account 
Coulomb interaction effects (in the sequential tunneling 
limit). Finally, in Sec. IV we give some final remarks. 
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FIG. 1: Sketch of the system: a quantum dot coupled to two 
ferromagnetic leads via tunnel barriers. The left FM lead has 
its magnetization fixed while the right-hand side can be either 
in parallel or antiparallel alignment. A pulsed bias voltage of 
duration s is applied across the system in order to generate 
transient spin-polarized currents. When the bias voltage is 
turned on (0 < t < s) the dot's level moves into resonance 
with the emitter states, and the dot becomes populated (a 
charging process) with a current passing through it. When 
the bias is turned off (t > s) ed moves above /i£ and [ir and 
the dot's occupation decays into the leads (a discharging pro- 
cess). Due to the ferromagnetism of the leads these transient 
charging and discharging processes become spin dependent. 



II. TRANSPORT FORMULATION 

To describe the system of a quantum dot coupled to 
two ferromagnetic leads, see Fig. [TJ we apply the follow- 
ing Hamiltonian 

\fL<JT) 

+Un 1 n l , (1) 

where eko-ij(i) is a time-dependent free-electron energy 
with wave vector k and spin a in lead r) (77 = L, R). This 
energy can also be written as e\s_ ari {t) = ^ct^ + ^W' w ^ n 
e k being the time-independent energy and A v (t) gives 
the time evolution of the external bias. The energy edit) 
is the time-dependent spin-degenerate dot level, which 
can also be written as ed(t) = e d + Ad(t), where e° is the 
time- independent level and Ad(t) follows the bias volt- 
age. It should be noted that in a quantitative theory one 
should consider a level-shift A^, which depends on the 
level occupation, via some suitable self-consistent proce- 
dure. We shall address this issue in our future work, but 
for the present purpose the simple model suffices. 



The operator Cko-ij (cj^^) is an annihilation (creation) 
operator for a single-particle momentum state k and spin 
(7 in lead 77 (77 = L,R), and d a {d) a ) is an annihilation 
(creation) operator for the single-particle dot's state Cd- 
The matrix element V\<_ ari ^ couples the leads with the 
dot, and we assume that the tunneling process is spin- 
independent. Finally, the [/-term describes the Coulomb 
repulsion in the dot, with n a = d\d a . 

In order to calculate the current we use the definition 
1% = —e(N%), where e is the electron charge (e > 0) and 

iVJ = J2k c t.a-q Ciicr v i s tne to t a l number of electrons with 
spin a in lead 77. From this definition it is straightforward 
to show tha1^^ 



12{t) = 2 e Re{^lWG< kCTi; (M)}, 



(2) 



where 



with GTjpfah) being the retarded (lesser) Green func- 
tion of the dot and frii^a) i s the time-independent Fermi 
distribution function of lead 77. Substituting Eq. © into 
Eq. © and following Ref. [38| we find 
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Q-iJt 1 dt2t kar ,(t 2 ) 
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(4) 



with T%(e,t lt t) = 2np an (e)\V an (e)\ 2 e l ^ dt2A ^ t2} . 
These results are exact, and they can in principle be 
used to study the intricate interplay between time- 
dependence, coherence and interactions. Their use, 
however, requires the knowledge of G r and G < , which 
come from the solution of the noncquilibrium Dyson 
and Kcldysh equations, respectively. For our main find- 
ings, though, it is sufficient to consider a non-interacting 
model, for which an exact solution can be obtained. 
Next, in Sec. 3D, we show that our results change only 
slightly when Coulomb interaction is included in a master 
equation based scheme. 

In the wideband limit (WBL) , 40 ' 41 and for noninter- 
acting electrons Eq. (H} can be written as 

JJ(t) = -eT*{Mt)) + J ^/^ImK.M)]}, (5) 

where (n„) is the time-dependent dot's occupation, given 
by 

(n a (t)) = lm{G< a (t,t)} 

= E r " / f^)IA, t ,M)| 2 , (6) 

and A ar/ (e,t) is defined as 



dtiG; CT (i,ti)e [ " (t -* l) - < ^ 1 dtA "©l. (7) 
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The retarded Green function in the nonintcracting model 
is given by 

G^(Mi) = -ie(t-*i)e-Tf(*- tl )e- <J "A« <Fe -© J (8) 

where r a = + T^. For a voltage pulse V(t) = 
Vo0(t)6(s — t) (see Fig. [IJ, and assuming that this pulse 
is applied on the right ferromagnetic lead, with a linear 
bias drop along the junction, we have Ai(t) = — Vl = 0. 
Afl(t) = -V R (t) = -V(t) and A d = -V d = -V(t)/2. 
With these definitions, we find for < t < s 4 ^ 



A aT] (e,0 <t<s) = 
- 

and for t > s we obtain 
A ar] {e,t > s) = 



e i{e-eo+V d -V v +iV,/2)t 

e-e Q + iT a /2 
2 _ e i(e-e +V d -V v +iT a /2)t 

e-e + V d -V ri + iT a /2 ! 



J(e-e +ir a /2)t i(V d -Vr,)s 



(9) 



at] y-i - - "j i -T-i /n 

e i(e-eo+ir„/2)(t-s) _ e i{V d -V v )s e i(e-e +ir „ /2)t 

e-e + V d -V v + iY a /2 

1 _ gi(e-e +ir CT /2)(t-s) 



e - e + ir CT /2 



(10) 



By substituting Eqs. © and ([13 into Eqs. ©-© yields 
the final result for the spin-resolved occupations and cur- 
rents. Numerical results are described in the next section. 



III. RESULTS 
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FIG. 2: Occupations n-f (solid line) and ri\ (dashed line) and 
the spin accumulation m = n-f — n\ (dotted line) as a func- 
tion of time for both (a) parallel and (b) antiparallel config- 
urations. When the bias is turned on (off) the dot is charged 
(discharged) in a spin-dependent manner. This results in a 
time-dependent spin accumulation, with a sign reversal in the 
parallel case. 



A. Parameters 

In our numerical calculations we assume that the volt- 
age pulse is applied to the right electrode, so that fiR = 
—V(t) while hi, is kept constant equal zero. The dot's 
level is taken originally (zero bias) above the chemical 
potentials \ii and hr> e o = 0-5 meV. The temperature is 
assumed to be T = 2.5K (fc^T « 215/ieV), thus allowing 
a small thermally excited occupation of the dot in equi- 
librium. To describe the ferromagnetism of the leads we 
choose the tunneling rates to be = F [l + (—l) Sl <*p] 
and r£ = r [l ± (-l) s i"p], where r is the leads-dot 
coupling strength and p gives the polarization degree 
of the leads4^ Here we assume a weak coupling with 
r = 1 ^teV ) 44 ' 45 and a polarization degree p = 0.4. The 
+ and — signs in give the parallel and antiparallel 
configurations, respectively. Due to the ferromagnetism 
of the leads (p ^ 0) we have Tf > Tf and Tf > Tf 
in the parallel case and the opposite T? < Tf in the 
antiparallel alignment. For the bias voltage we adopt 
V(t) = V 8(t)8(s - t) where V = 5 mcV and s = 3 
ns4£ The charging energy U is set equal to zero in Sec. 
III(b)-(c) and equal to 3 meV in Sec. 111(d). 



B. Spin-polarized occupations 

Figure [5] shows the spin-resolved occupations nj and 
nj_ and the spin accumulation m = nj — as a function 
of time for both (a) parallel and (b) antiparallel configu- 
rations. Before the bias is turned on the level e d is above 
the electrochemical potentials \if\ ( ? 7 — L,R), and the 
dot is only slightly occupied due to thermal excitation. 
When the bias is turned on at t — the dot's level is 
brought into resonance (/zj, < e d < /j,r), thus resulting 
in an enhancement of n a and m. In the parallel case 
[Fig. EJa)] the spin up population increases faster than 
the spin down one, and both attain the same stationary 
value around 0.5. The steeper enhancement of n\ com- 
pared to n\ is related to the inequality Ti 1 > Tf, that 
gives a faster response for the spin j component. Since 
= in the P case, the in- and out-tunnel rates 
compensate each other, thus resulting in n-f = for 
asymptotic times. When the bias voltage is turned off, 
e d raises above \il and [ir and the population of the dot 
begins to decay, with a faster discharge for the \ com- 
ponent. The spin accumulation reflects the dynamics of 
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n-j- and rij. In the range < t < s, m reaches a local 
maximum due to the faster enhancement of compared 
to 7ij_. In contrast, when the bias voltage is turned off 
(t > s), m shows a local (negative) minimum due to the 
fast discharge of . 
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FIG. 3: Spin resolved currents against time for both left and 
right leads and both alignments. In both configurations the 
currents in the left (right) lead are suppressed (enhanced) just 
after the voltage is turned on, and then they attain stationary 
plateaus. When the bias voltage is turned off (t > 3ns) the left 
and right currents become the same for each spin component 
in the P configuration, while in the AP case the f current 
becomes bigger than the J, current in the left lead. 



In Figure[5Jb) we show the evolution of the occupations 
and the spin accumulation in the antiparallcl alignment. 
We note that rif increases faster than nj, as in the P case. 
In contrast, though, n.j attains a higher value than n± in 
the stationary regime. This is related to the out-tunnel 
rates which are now inverted with respect to the parallel 
case: Tf < Tf. When the bias is turned off both n-f and 
n-i decrease due to the transient discharge. In particular 
the spin up electron population discharges predominantly 
to the left lead while the spin down component discharges 
to the right, following their corresponding majority den- 
sity of states (or equivalently the majority tunnel rates). 
The way how spins f and j. charge and discharge are more 
clearly seen in the spin-resolved current curves described 
in the next section. 
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FIG. 4: Total current in the left ferromagnetic lead, If + 



7i , for both parallel (solid line) and antiparallel (dotted line) 
configurations. After the bias voltage is turned off (t > 3 
ns) the total antiparallel current becomes greater than the 
parallel one, lifted by the spike of j current seen in Fig. IS^b) . 
This results in the time-dependent negative TMR seen in the 
inset. 



C. Spin- resolved currents 

Figure [3] shows if and I± for both leads and both ferro- 
magnetic alignments. In the P configuration [Fig. [3ja)] 
the left currents If and If show a transient suppres- 
sion and then attain their respective stationary values 
with if > if. In the right lead the currents if and if 
increase (in modulus) up to their respective stationary 
values. When the bias voltage is turned off 1% becomes 
negative as iff. The negative sign of both 1% and I„ 
means that the electrons are flowing from the dot to the 
leads (discharge). In particular the spin J, electrons dis- 

L R L R 

charge much slower than the f ones, due to T ^ ' < ' . 

In the AP configuration [Fig. \3[h)] if and if show 
a suppression just after the bias voltage is turned on, 



then they attain a stationary value with if 



If. In the 



right lead the currents if and If are enhanced until they 
reach equal plateaus. When the bias voltage is turned 
off, if and if change sign (discharge of the dot) and a 
spike of spin f current in seen in the left lead (if ^ if). 
This reflects the preferential discharge of spin up elec- 
trons to the left lead, according to Tf > Tf. No spike 
is seen in the parallel configuration, where spin up elec- 
trons discharge equally to both leads. In contrast, in the 
AP alignment the spin down electrons discharge prefer- 
entially to the right lead due to the inverted inequality 
rf < Tf, while in the P case its discharge is equally to 
both sides (Tf = rf). 

Negative TMR. In figure fj| we show the total current 
if) for both parallel and antipar- 



in the left lead (If 
allcl configurations. 



Due to the strong spin-polarized 
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discharge (t > 3 ns) in the left lead when the sys- 
tem is AP aligned, the total current obeys the unusual 
inequality l\ P > Ip, which results in time- dependent 
negative tunnel magnetoresistance (see inset), defined as 
TMR = (J| - Iap)/Iap- As the time evolves the TMR 
keeps increasing, due to the longer spin-down lifetimes 
when the system is parallel aligned. More specifically, in 
the AP configuration both spin up and down discharge 
fast to the left and to the right leads, respectively, fol- 
lowing their majority spin populations (or equivalcntly 
the tunneling rates). In contrast, in the P alignment 
the majority populations occur for spin up in both leads 
(T^' R > T^' R ). This turns into a fast discharge for spin 
up electrons and a slow discharge for the down compo- 
nent. This slow spin down discharge sustains the total 
current much longer than in the AP configuration, and 
eventually for long enough times we find Ip 3> l\ P - 

Displacement Current. In the transient regime the left 
and the right currents are not in general the same (II ^ 
/r), due to charge accumulation/depletion in the dot. 
The generalized conservation law is given by the continu- 
ity equation I^+I R — l£ ls = 0, where is the displace- 
ment current for spin a, given by I% ls = ed(n a (t)) jdt. In 
order to check the accuracy of our numerical calculation 
we have verified numerically the continuity equation. 



D. Effects of Coulomb Interaction 

An exact treatment of the Coulomb interaction rep- 
resents a formidable problem, and in the context of the 
present Hamiltonian only few results arc known in equi- 
librium, and none in nonequilibrium, even less so under 
transient conditions. Nevertheless, in certain limits ap- 
proximate treatments may give a good qualitative un- 
derstanding of the generic behavior. One such case is the 
sequential-tunneling limit (ro <C ksT), where the Master 
Equation (ME) approach is known to work well. Here, 
we use the ME to estimate the effects of Coulomb inter- 
action in our results4i The current expression is given 
by^ 



(b)-antiparallel 



(1-/ 1) )P 2 ], (11) 

(n CT (l - ng-)) and 



where P = ((1 - n T )(l - n^)), P a 
P2 = {n-\ni), are the probabilities to have no electron, 
one electron with spin a and two electrons, respectively. 
The Fermi functions f n and are evaluated at e c i and 
f-d + U , respectively For the dot's occupation we write 



dt 



1 



it 



= Y^muPo - (1 - f v )p* 



+f v P*-(l-f n )P 2 ], (12) 
and for the double occupancy probability we have 

j t {^H) = £r2[/„p ff - (1 - f v )p 2 }. (13) 




err) 



2 3 4 5 

time (ns) 



FIG. 5: Spin-resolved occupations (a)-(b) and currents (c)- 
(d) in both parallel (P) and antiparallel (AP) configurations. 
We take (7=1 meV and U — 3 meV, labeling the traces 
by (1) and (3), respectively. For U — 1 meV the results 
are indistinguishable from the U — case. For U — 3 meV, 
though, we find a suppression of the spin-resolved occupations 
and currents. In the AP alignment this suppression results 
in an enhancement of the spin accumulation and in a spin 
polarized current in the stationary plateaus (|7^ ,i? | > 
To clarify the range t > 3 ns, in panel (c) the currents 7^ 
are on top of 1^ while in panel (d) they are apart from each 
other. In addition in the AP case, \I^\ and for U = 1 
meV are slightly greater than their corresponding values for 
U = 3 meV, and 7^ is almost on top of 7f*. 



For the nonintcracting case (U = 0) the time- 
dependent results obtained from Eq. (jTTJ) are identical 
to those seen in Sec. III(b)-(c). For the interacting case 
(U 7^ 0), we find that for U =1 meV the results are 
indistinguishable from the U = case [see Fig. [5]. This 
is so because for small enough U both channels td and 
£d + U attain resonance for V(t) = 5 meV. In contrast, 
for U — 3 meV the channel e<j + U remains above the 
emitter chemical potential when the bias voltage is ap- 
plied, which turns into a suppression of the occupations 
and the currents. In particular in the AP configuration 
this suppression is stronger upon the spin down compo- 
nent, seen in both occupations [panel (b)] and currents 
[panel (d)]. This is due to the spin imbalance n-\ > ni 
typically present in the antiparallel alignment. This spin- 
polarized suppression in the AP configuration gives rise 
to an enhancement of the spin imbalance [see Fig. E)(b)] 
and to a spin polarized current (|7^ ,i? | > |) in the 
stationary plateau. 

In Fig. ^ we see the effects of U on the dynamical 
TMR. For U = 1 meV the TMR is basically the same 
as before [Fig. gfinset)]. For U = 3 meV the TMR 
is enhanced (in modulus) for both on and off voltage 
regimes (0 < t < 3 ns and t > 3 ns, respectively). In 
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FIG. 6: Tunnel magnetoresistance (TMR) as a function of 
time for U = 1 and U = 3 meV. For U = 3 meV the TMR is 
enhanced in the stationary regime (2< t <3 ns) and becomes 
even more negative after the bias voltage is turned off (t > 3 
ns). 



particular the Coulomb interaction turns the TMR even 



more negative after the bias voltage is turned off, which 
reaches —40 % around 3.5 ns for U = 3 meV. 



IV. CONCLUSION 

We predict novel spin-dependent effects in a quantum 
dot coupled to two ferromagnetic leads driven by a rect- 
angular bias voltage pulse. Based on nonequilibrium 
Green function and master equation techniques we cal- 
culated the spin-resolved occupations and currents, the 
spin accumulation and the tunnel magnetoresistance in 
the transient just after the bias voltage is turned on and 
off. Our main findings are: (i) a sign change of the spin 
accumulation as the time evolves in the P configuration, 
(ii) a spike of spin f current in the emitter lead when the 
system is antiparallel aligned, and (iii) a time-dependent 
TMR that attains negative values. This negative amount 
can be further enhanced due to intradot Coulomb inter- 
action. 
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